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Abstract. We study the interface between two different Z{3) vacua in the decon- 
' fined phase of SU (3) pure gauge theory in 2 + 1 dimensions just above the critical 

temperature. In simulations of the Euclidean lattice gauge theory formulation of the 
^ . system we measure the fiuctuations of the interface as the critical temperature is 

approached and as a function of system size. We show that the intrinsic width of the 
interface remains small even very close to the critical temperature. Some dynam- 
ical exponents which govern the interaction of the interface with our Monte Carlo 
algorithm are also estimated. We conclude that the Z{3) interface has properties 
broadly similar to those in many other comparable statistical mechanical systems. 
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1 Introduction 



It has been known for some time that pure gauge theories have a non-confining 
high-temperature phase [I[]and so any theory which is in a confining phase at zero 
temperature (T = 0) must have a deconfining phase transition at some finite critical 
temperature, Tc. Pure non-abehan SU{N) gauge theories, such as QCD without 
dynamical fermions, display this behaviour in both 2 + 1 and 3 + 1 dimensions. At 
low temperatures the colour charge of QCD is confined which is why no free quarks 
are seen in the relatively cool universe of today. At high temperatures, there is a 
non-confining phase, corresponding to the free quark-gluon plasma believed to exist 
in the hot early universe. In fact there are in general several degenerate phases of 
this type - as many as the number of elements of the centre, Z{N), of the gauge 
group. Fermionic matter breaks the vacuum degeneracy, albeit on a small scale, so 
we consider only the pure gauge theory in what follows. 

Two different types of interface are possible in the theory. The first is between 
the ordered and disordered phases and is only stable at the critical temperature 
where the two phases can coexist. Secondly, an interface can form between two of 
the ordered phases with different Z{N) vacua, and this is usually called an "order- 
order interface", or ^^Z{N) interface". Clearly, this type of interface can only exist 
above the critical temperature. 

From the thermodynamic point of view there is something strange about these in- 
terfaces. On entropic grounds one expects to find disorder at high temperatures and 
order at low temperatures, not the other way round. The fundamentally Euclidean 
nature of this picture may be responsible |^, §|; there is no sensible counterpart 
to the Euclidean order parameter in Minkowski space (a naive Wick rotation leads 
to an imaginary time-like gauge field in Minkowski space) and so it may be that the 
interfaces do not exists as physical objects in the high temperature universe. In- 
deed, it has even been suggested that only one true physical phase exists, even in 
Euclidean space, at high temperatures but recent numerical work does not support 
this point of view. Whether or not they are a Euclidean artifact, by virtue of their 
contribution to the partition function and to expectation values calculated using the 
Euclidean path integral, the interfaces must be included (like the instanton) in a 
non-perturbative analysis of the thermodynamics of QCD. 

Recently the properties of Z{N) interfaces at very high temperatures have been 
investigated numerically in 2 + 1 dimensions 0, H. These simulations give a non- 
perturbative check on the conclusions of analytic calculations [1^ and show very 
good agreement; thus there is no reason to suppose that the Z{N) phases are not 
distinct at high temperatures. In this paper we investigate the behaviour of the 
interface at much lower temperatures close to the critical temperature where it dis- 
appears. We shall be concerned not with the usual thermodynamical quantities. 



which have been thoroughly investigated |11|, but rather with the geometrical prop- 



erties and structure of the interface itself. 

This paper is organised as follows. In section two we review briefiy the theory of 
the deconfinement transition and the theoretical basis on which we shall analyse the 
simulation data on the interfaces. In section three we discuss the practical difficulties 
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in identifying the interface from the gauge field configuration and describe how we 
have overcome them; some other aspects of the data analysis are also discussed. 
Section four contains the results that we have obtained characterising the interface 
fluctuations, and section five gives our results for various other properties of the 
interface. In section six, we give our conclusions and discuss some of the remaining 
puzzles associated with these interfaces. 



2 Theoretical Background 

In the Euclidean formalism, equilibrium finite-temperature field theory is obtained 
by considering a system which is of infinite volume but compact in the Euclidean 
time direction with period /3t = T~^. In a lattice field theory implementation, it is 
of course not possible to vary the number of lattice spacings in the time direction, 
Lt, continuously; instead we work with a fixed Lf and vary the gauge coupling, /3, 
so that the lattice correlation length, and hence the physical value of the lattice 
spacing, a, and hence the physical Pt = Lta, vary continuously. 

The phases of a gauge theory at finite temperature are characterised by the free 



energies of static configurations of quarks and antiquarks [1^. The self-energy of 
a single quark in the gluon medium, Fg, is related to the expectation value, L(x.), 
of a single Polyakov line, a time-ordered Wilson loop which wraps around the time 
boundary for a fixed spatial location x. 



(x) = iTrre^^/o^^'^"^o("'^) = e-f^^^" 



If < L{x.) >= then the insertion of a single quark will require infinite energy, 
corresponding to a confining phase. In contrast, if < L{x.) >^ then the insertion 
energy will be finite and isolated colour charges can exist - this is the deconfined 
phase. Thus, the expectation value of the Polyakov line gives an effective test for 
confinement. 

A general non-Abelian gauge transformation which leaves the Euclidean QCD 
action, Se, invariant will also leave L invariant if the transformation is periodic 
in time. However, Se is actually invariant under an additional global symmetry. 
In the lattice formulation this consists of multiplying each time-like link variable 
at a given time slice by an element of the centre of the gauge group (the set of 
elements which commute with all members of the group, Z{N) in the case of SU (N)). 
This global symmetry is an invariance of Se and cannot be undone by a gauge 
transformation. The topologically non-trivial Polyakov line, which wraps around 
the periodic boundary condition in the time direction, is not invariant under these 
transformations, but is rotated by an element, z, of the centre: -^^(x) — » zL{x.). 
Clearly, < -L(x) > acts as an order parameter for the centre symmetry, distinguishing 
between broken and unbroken phases, and between the different broken phases 

g. 

In the deconfined phase, where the Z{N) symmetry has been spontaneously 
broken, the system can exist in any one of degenerate vacua. It is possible to 
arrange boundary conditions of the system so that different parts of it exist in 
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different vacua, i.e. different Z{N) phases. The existence of these distinct domains 
forces the appearance of domain walls, or ^^Z{N) interfaces", where they meet |P, 10|. 



Within these interfaces, the gauge fields interpolate between the different vacua, as 
does the expectation value of the Polyakov line. A perturbative calculation at large f3 
leads to the conclusion that an instanton interpolating between two vacua should 
have characteristic size, Iw of order l/gy/T; therefore, at least at large enough /3, 
we expect the interface to have an intrinsic width of this magnitude as well. This 
length scale is the same order as the Debye screening length. Id which governs 
gauge-invariant correlation functions of the time-like component of the gluon field 
{Aq) at large distances and high temperatures. The free energy of a quark- ant iquark 
pair, over and above the sum of their separate free energies, vanishes as the quarks 
become infinitely far apart, and, from perturbation theory, we know that quarks are 
screened at a distance of the order of the inverse electric mass, so rather than being 
logarithmic in x, the interaction potential takes the form 



F,g(x)-2Fg = l^(|x|,T) ~ -Cexp -2^ for T > T,. (2) 

|x|^oo y L]J J 

The factor of two arises because gauge invariance leads to an exchange of two gluons 
being the lowest-order contribution in perturbation theory. The Debye screening 



length is known |]14| to have logarithmic corrections so that 



'°-9VT-4log(i|) + ... 

The corresponding corrections to Iw have not been calculated so, even at very high 
temperature, a direct comparison between the interface width and the Debye screen- 
ing length is not possible. At low temperatures, near the critical point, we expect all 
these quantities to be functions oi (3 — (3c and there are no analytic results available 
with which to compare our simulations. 

At very high temperatures the magnitude of < L > is large; the change in its 
argument from one Z{N) phase to the other is abrupt. The interface between the 
two phases is sharp (see fig. 2 of section 3); it is essentially a one-dimensional object 
that separates the phases and exhibits relatively small transverse fluctuations. As 
(3 decreases towards its critical value, < L > vanishes according to 

<L>~(/5-/5,)^-, (4) 

where /3c = 8.175(2) in the infinite spatial volume limit on A^^^ = 2 lattices, and the 
measured critical exponent Pm is consistent with the universality prediction 
of I from the three-state Potts model [|l^. However, this only tells us about the 
expected height of the interface and not about its geometrical nature as /? decreases 
(except that it must ultimately disappear at f3c). It may be that, as the free energy 
per unit length declines, the interface undergoes increasingly violent fluctuations 
while remaining intrinsically a one-dimensional object, the height of the interface 
declining until it vanishes at jSc. An alternative is that the interface may become 
intrinsically broader, thus becoming more two-dimensional, with more and more 
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Figure 1: The vacuum diagrams corresponding to the 2nd, 4th and 6th connected 
vacuum correlation functions respectively. The open circle simply denotes point x, 
the endpoint of each external leg; the solid circle denotes an interaction vertex, A; 
and the momentum travelling along a leg is denoted by k. 

regions of disordered phase within it. At f3c these disordered regions expand to take 
over the whole system and the interface disappears. 

To help in the analysis of the fluctuations of the interface we shall compare them 
to a simple model. Our lattice is taken to have ~ 3Ly ^ ; the boundary 
conditions are periodic and a "twist" is introduced (after to ensure that for 
X ~ the system is in one Z{3) phase while for x ~ Lj, it is in a different one. The 
interface thus forms stretching across the lattice in the y direction and separating 
the two phases. Now let (j)(y) denote the displacement, in the x direction, of the 
interface from its equilibrium position and construct a simple effective Lagrangian, 
C, for (f). Translational invariance dictates that C can only depend on and 
invariance under — >■ —(j) rules out odd powers; if the interface were just executing 
a free random walk across the lattice then C ~ ((9^0)^. To obtain non-gaussian 
moments a correction term must be added and the simplest resulting Lagrangian is 

which fits the data quite well with 

A = iXo 

Ao ~ {f3-f3cr (6) 

where the length scale i appears in A on dimensional grounds. It can be any combi- 
nation of the two length scales which may be relevant, Ly and Id- The diagrams for 
the first three even moments of are shown in fig.l. It is straightforward to com- 
pute the connected correlation functions from C although the mode sums have to 
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Figure 2: Polyakov line (real part) profiles of the same Z(3) interface at different 
temperatures. The top left picture shows the interface on a 2 x 24 x 96 lattice 
with f3 = 50. Going from left to right, and then top to bottom, the other pictures 
illustrate an interface between the same phases (< L >^ ^^m/s ^j^g |g£^ ^j^g 
< L >?» 1 on the right) for /3 = 10, 9, 8.5, 8.3 and 8.1 respectively. 

be done numerically. Putting £ = Ly the connected correlation functions (denoted 
by < . >c) are then predicted to behave like 

<ct>'>^ ^ Lir'xl (7) 
3 Identifying the Interface 

Fig. 2 shows the ^(3) interface on a particular gauge field configuration for a number 
of (3 values. As discussed in section 2, at high (3 the interface is sharp and well 
defined; however by /3 = 9, still some way above the critical value, it is becoming 
difficult to identify the location of the interface by eye with any degree of certainty. 

One way to deal with the increase in phase fluctuations is to cut out the highest 
frequency modes altogether. This can be achieved most simply by a "box-car" 
average, where the Polyakov lines of each three-by-three spatial array of points are 
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Figure 3: Smoothed Polyakov line profiles, produced by processing the raw pictures 
of fig. 2. As in that figure, P — 50, 10, 9, 8.5, 8.3 and 8.1 respectively. 

averaged over to give a new value, allocated to the central point. The effect of 
smoothing out the fluctuations is shown in fig. 3, produced by applying the box-car 
average to the pictures of flg.2. To make things even clearer, the Polyakov lines 
have first been re-binned so that the 24 x 96 spatial dimensions become 12 x 48, 
with each two-by-two array being averaged to one point. Now the Z{3) interface 
can clearly be seen for all but the last picture, i.e. for all temperatures above (3c. 
The fact that the interface can still be seen for temperatures just above the critical 
value suggests that a study of its structure should still be possible quite close to the 
phase transition. 

However, there is a drawback with this smoothing technique as it averages out a 
whole class of fiuctuations which may be important in the interface collapse. A more 
satisfactory method is to produce a contour map from the real part of the Polyakov 
lines on the lattice. Consider following a contour whose height is mid-way between 
that of the average Polyakov line at the low- a; end of the lattice and that at the high- 
X end. This picks out the mid-height of the interface as it goes across the lattice, as 
well as any bubbles of fiuctuating phase which are large enough to cross the contour 
height. The crucial point to note is that only the Z{3) interface will cross the entire 
lattice, i. e. the only contour that will wrap once around the lattice in the transverse 
(y) direction is that corresponding to the interface. Any contour representing a 
bubble of phase will join up with itself without a net crossing of any lattice boundary. 
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(i) 




(ii) (iii) 

+ + + 
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Figure 4: Options for the contour-following algorithm. Given a point at the desired 
contour level on the lower side, the Polyakov lines at the upper two corners determine 
the progress of the contour in cases (i)-(iii). A plus (minus) sign indicates a Polyakov 
line value greater (less) than the contour level. The same corner configurations of (iv) 
and (v) give two possible contour directions; we choose between them by considering 
the central average of the corner values. These pictures exhaust all possibilities for 
the values of the lower corners given, but any general configuration can be rotated 
and reflected into one of these. 



This enables us to tell the interface apart from any fluctuations, and to identify the 
position of its mid-height. By following contours at various heights between the two 
extremes, we can study more of the structure of the interface. The advantage of 
using contours is that no prior processing of the Polyakov line data is required. We 
can also study contours of the box-car- averaged data and this can be useful very 
close to the critical temperature. 

To flnd a contour at a particular level, we return to the raw Polyakov line data. 
For y = 0, the algorithm follows a procedure similar to that above: starting from 
the centre, x — Lx/2, it looks for the nearest point passing through the contour 
level. Then, it considers the square of points formed by the two neighbouring lattice 
sites at y = and the corresponding sites at y = 1. The algorithm decides to which 
side of the square the contour should go, as detailed in fig.4. Then, starting with 
this side, the whole procedure continues. Eventually, the contour must join up with 
itself. Repeated application of this technique for different contour levels determines 
the structure of the interface to the desired precision. The same procedure can be 
applied to the smoothed Polyakov line conflgurations, to obtain smoothed contour 
maps. An example of a contour map is shown in fig. 5, for the same /5 = 10 gauge 
configuration featured in fig. 2 and fig. 3. Note that the interface in the centre of 



8 



Contours of Raw Polyakov Lines 



c 

o 



> 




40 60 
Longitudinal Direction, x 

Contours of Smoothed Polyakov Lines 




40 60 
Longitudinal Direction, x 



Figure 5: Contour maps of raw and processed Polyakov lines for (3 = 10 on a 
2 X 24 X 96 lattice. These maps correspond to the (3 = 10 pictures of fig. 2 and fig. 3 
respectively. The contour levels are at —0.25, 0, 0.25, 0.5 and 0.75, darker regions 
corresponding to higher values. 
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the maps is represented by the only contours going all the way across the lattice. 
The first picture shows the raw data of fig.2; the second, the smoothed data of fig. 3. 

Left to its own devices, the interface will execute a random walk in the x direc- 
tion. When the interface reaches x ~ or x ~ it will interact with the twist 
and tunnel into different Z{3) vacua, producing a different Z{3) interface each time. 
However to study the profile of one particular Z{3) interface without interference 
from other vacua the interface should be held in place well away from the twist, 
so that no part of it, even its wildest fluctuation, passes through the twist. After 
each Monte-Carlo (MC) sweep of the lattice, we locate the position of the interface, 
and then slide the whole gauge configuration along the lattice until the interface is 
re-centred at x ~ Lx/2. Any link variables which are shifted through the twist are 
multiplied by the appropriate factor, so there is no physical effect and no violation 
of translation invariance arising from this procedure. 

For our MC simulations, we used a mix of Cabibbo-Marinari heat-bath steps 



Tql , using the Kennedy-Pendleton algorithm |jT9| to update SU (2) subgroups, and 



Brown- Woch over-relaxation steps after using many initial heat-bath sweeps to 
equilibrate the system. Measurements of physical operators took place only every 
four sweeps, to reduce the correlations caused by adjacent configurations. Since 
we were interested in the evolution of the interface, contour measurements were 
taken after every sweep of the lattice. In order to study finite-length effects for the 
interface, we performed simulations on a number of different lattice widths [Ly = 
18, 24, 30, 36, 42, 48, 54) and length = 72, chosen after studies with various lengths 
to ensure that not even the wildest fluctuations of the interface would reach as far as 
the twist during our full-length runs. We used (3 values of 8.25, 8.35, 8.50, 8.75 and 
9.00 to study the behaviour of the interface close to (3c- The simulations consisted, in 
each case, of 2k heat-bath sweeps and 100k main sweeps. Approximately 100 — 200 
hours of CPU time on a DEC 2100 A500MP machine were needed for each value of 
(3. 

The statistical analysis of our data, and especially of the Wick-subtracted mo- 
ments which are not expected to be normally distributed, is quite complicated and 
requires a method which does not rely on a priori knowledge of their distribution. 
We used methods based on the bootstrap principle and full details are given in 
l2l. 



4 The Fluctuation Moments of the Interface 

In fig. 6, a contour snapshot of the interface is shown for each temperature studied, 
at the same point in each simulation, just over halfway through. Three contours 
have been followed in each case. We can see that the contours stay close together 
across most of the interface, indicating that the interface remains relatively narrow 
as the temperature drops. The shape of the interface appears to fluctuate rather 
more wildly with the falling temperature, but the interface does not spread out 
noticeably in width. 

As discussed in section 2, the fluctuations are characterised by their moments 
< >c] 4'{y) can be the displacement from its mean position of any contour 
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P=8.35 
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P=9 Contour levels: -0.059, 0.119, 0.296 



P=8.75 Contour levels: -0.056, 0.112, 0.279 



P=8.5 Contour levels: -0.051, 0.102, 0.256 



P=8.35 Contour levels: -0.042, 0.097, 0.235 



P=8.25 Contour levels: -0.037, 0.088, 0.214 



Figure 6: Contour maps of the interface on a 2 x 36 x 72 lattice, for decreasing 
temperature. Moving left to right, and then top to bottom, the picture illustrates 
the interface structure around the longitudinal centre of the lattice after 2k heat- 
bath and 60k main sweeps for /? = 9, 8.75, 8.5, 8.35 and 8.25. 
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that has been identified. We measured the first six moments, and took separate sets 
of data for each of the three contours followed, at 25%, 50% and 75% of the height 
of the interface respectively. The qualitative behaviour is roughly the same for each 
contour. The upper and lower contours are subject to rather more fiuctuation than 
the middle one, as we expect since they are closer to the outside of the interface 
and therefore more susceptible to bubbles of phase forming nearby and distorting 
the edge of the interface. The even moments behave very smoothly, whilst the odd 
moments are very much smaller and show more variation; the first moment must 
average to zero by definition, and we expect all odd moments to tend to zero for 
infinite statistics. We also computed the behaviour of the moments of smoothed 
data and found that it is very similar to that of the raw data, both qualitatively 
and quantitatively. This indicates that the smoothing procedure does not lose a 
significant amount of information about the fluctuations and is consistent with the 
notion that the important fluctuations take place on a distance scale much larger 
than the lattice spacing. 

Our results for jS > 8.35 show reasonable statistical distributions for the moments 
over the 100k sweeps. However, those for f3 = 8.25 and smaller Ly occasionally con- 
tain "spikes", observations for the moments which are far larger than we would 
expect, sometimes by several magnitudes. The spikes are not caused by wild fluc- 
tuations passing through the twist, since similar spikes are not seen for the widest 
lattices where fluctuations are larger. They are an artefact of the greater energetic 
instability of the smaller interfaces near the critical point. If parts of the inter- 
face collapse during interaction with nearby bubbles of phase, the interface contour 
may briefly undergo a drastic distortion or partial collapse, even appearing to pass 
through the twist because of the change in Polyakov values there and confusing our 
algorithm. We cannot prevent the spikes occurring, since the cause is a real physical 
process, but they can have a serious effect on the overall averages, especially for 
the higher moments. Thus, we would like to be able to examine the data with the 
spikes removed. To do this, we impose an upper cut-off on the statistical distribu- 
tions. To ensure that only the spikes are removed, we need to deflne the cut-off in 
a statistically sensible manner. We take it, for each distribution, to be the median 
plus 2.5 times the interquartile range of the distribution of logarithmic values of 
data, where the interquartile range is defined to run between the 25th and 75th 
percentiles of the distribution. This ensures that the cut-off corresponds only to 
the upper tip (to be precise, the top 0.04%) of a normal distribution, ensuring that 
only data sets containing spikes are affected by the procedure. The exact multiplier 
"2.5" can, of course, be varied, but we find that a lower number {e.g. 2.0, cutting 
off 0.3% of the distribution), has too great an effect on the normal contribution to 
the distributions, whereas a higher number {e.g. 3.0, cutting off only 0.002% of the 
normal distribution) permits too much influence from the spikes. The procedure is 
fully detailed in p2|. In our tabulated results values which were obtained with this 



procedure in effect are marked with a f. 
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Table 1: Estimates of < 0^ > from Raw Data (Corrected), 50% Contour 
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Table 2: Estimates of < (j)^ >c from Raw Data (Corrected), 50% Contour 
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Table 3: Estimates of < 0^ >c from Raw Data (Corrected), 50% Contour 
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Figure 7: The connected moments a.i (3 = 8.5 plotted against the lattice width. 
Diamonds and boxes represent raw and smoothed data respectively. The lines are 
best fits to (8). 

Fig. 7 shows the dependence of the moments, calculated for both the raw and the 
smoothed 50% contour, on the lattice width at /3 = 8.5. The slopes for the second 
and fourth moments are in reasonable agreement with the simple model results 
although the agreement is not good for the sixth. In fig. 8 we show the same data 
but with spike suppression in operation. The results are very similar to those shown 
in fig.7 demonstrating that the spikes are not a significant problem at /3 = 8.5. The 
roughness exponents defined by 

are tabulated in table 4 for all the j3 values we have studied (an asterisk means that 
no statistically significant determination was possible). 

In fig. 9, we plot the moments at fixed Ly = 36 against (3 when the spike- 
suppression procedure discussed above is not operating. The lines are fits of the 
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Figure 8: The connected moments with spike suppression at (3 — 8.5. 



Table 4 



p 


2a2 


4^4 




8.25 


1.33(2)t 


2.12(7)t 


1.46(22)t 


8.35 


1.27(l)t 


2.17(6)t 


2.38(29)t 


8.5 


1.15(1) 


1.95(5) 


1.82(26) 


8.75 


1.05(1) 


1.54(4) 


* 


9.0 


1.03(1) 


1.54(3) 


* 
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Figure 9: The connected moments on the Ly — 36 lattice plotted against p. The 
lines are best fits to (9). 
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Figure 10: The connected moments with spike suppression on the Ly = 36 lattice 
plotted against jS. The lines are best fits to (9). 

form 

<r>c--if3-f3c)-^" (9) 
taking (3^ = 8.175. Comparing the second and fourth moments with we find that 

^ ^ i/3 - /3^r-''('\ Ao ~ (/5 - (10) 

which then predict that 

<0^>,~(/3-/3J-^-^^(^i), (11) 

within roughly one standard deviation of the observed behaviour. In fig. 10 we show 
the same data as for fig. 9 but with spike suppression. We see that the results for the 
second moment are scarcely affected and the change for the fourth moment is fairly 
small but the effect on the sixth moment, which we would expect to be most affected 
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Table 5 



Ly 


72 


74 


76 


18 


.41(1) 


.99(6)t 


1.98(ll)t 


24 


.37(2) 


1.10(7)t 


2.10(15)t 


30 


.44(2) 


1.28(7)t 


2.58(25)t 


36 


.43(3) 


1.25(ll)t 


1.96(60)t 


42 


.44(2) 


1.29(5) 


2.09(12)t 


48 


.45(3) 


1.26(12) 


2.79(80)t 


54 


.46(4) 


1.32(11) 


1.87(57)t 



by the spikes, is substantial. Again comparing the second and fourth moments with 
(|^ we find that 

7~(/?-/3.)°-^^(^), A~(/?-/?e)°-^^(^^). (12) 
which in turn predict that 

<0«>c~(/5-/3e)-^-°«(^^\ (13) 

extremely close to the observed behaviour. The other sets of revised data now give 
very consistent results for 7 and A. 

Our other sets of data give similar results and the 7^ are tabulated in table 5 for 
all the L„ values we have studied. 



5 The Interface Width, Screening and Mobihty 

We define the transverse averaged Polyakov line by 

L{x,Ly) = L-^^Y.L{x,y) (14) 
y 

At very high temperatures where the interface is essentially rigid we expect it to 
show the classic instanton shape and be independent of Ly 

L(a;,L^) =/(x) ~ A tanh(^^^) +5 (15) 

where A and B are constants and / is the characteristic width of the interface. 
As the temperature decreases and the interface fiuctuates more (roughens) some of 
the spread in L{x) will be caused by the roughening; provided the fiuctuations are 
not too great, and assuming that the moments of are given by the free random 
walk results (good enough at the accuracy of measurement that we have), we can 
approximate the measured L{x,Ly) by 

L{x,Ly) = ^ ^ < /(x - 0(1/)) 
y 

= f{x) + Lyh{x) + Llh{x) + ... (16) 



18 



0.6 r 



0.4 



A 
J 

V 



0.2 



0.0 



-0.2 



-0.4 L 




25 30 35 40 45 

X 




Figure 11: The transverse Polyakov lines at = 9.0 and the weighted first differences 
discussed in the text. 
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Figure 12: The interface width deduced from transverse Polyakov lines plotted 
against 
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Interface Widths 




Figure 13: The average interface width, defined as twice the average separation of 
the 25% and 50% contours, measured in units of the lattice spacing. 

The piece hnear in Ly can be removed by taking weighted first differences of mea- 
surements at different Ly, 

L2L{x, Li) - LiL{x, L2) 



L2 — Li 

f{x) + 0{Ll,) (17) 



In the same way both hnear and quadratic Ly dependence can be removed by taking 
weighted second differences. Fig.ll shows the results of this procedure at /? = 9.0; 
once the part linear in Ly is removed all the curves collapse quite convincingly onto 
one. Fitting (|I^) to this universal curve a value for / can be extracted. As [3 is 
decreased the procedure, which essentially treats the fluctuations as a perturbation, 
works less well but the behaviour is reasonable down to (3 = 8.35; it was not possible 
to measure / at /3 = 8.25. The results for I are plotted in flg.l2 and show a mild 
increase as (3 decreases but we do not have enough information to extract any critical 
behaviour. 

We have also measured the interface width by using the contour data directly. 
The results of deflning the width as twice the separation between 25% and 50% are 
shown in flg.l3. Unfortunately we are only able to measure the separation in the 
X direction; as the interface fluctuates this quantity is bound to be greater than 
the shortest distance between the contours and so it is inevitable that, especially 
at the smaller j3 values, it will be over-estimated. Even without taking this into 
account we see that the width does not grow rapidly as the temperature decreases 
towards its critical value and the two methods for measuring the interface width 
yield qualitatively consistent results. 
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Figure 14: The wandering of the interface, for temperatures near the critical tem- 
perature, over 40k sweeps after an initial 400 heat-bath sweeps. 
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Figure 15: The divergence of the speed of the random walk in fig. 12 as the temper- 
ature approaches its critical value. 
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Finally, as a curiosity, we have examined some quantities which describe the 
dynamical (in computer time) behaviour of the interface and therefore relate to the 
efficiency of our MC algorithm in equilibrating the interface itself. The speed of the 
random walk in the x-direction made by the interface can be measured by turning off 
the mechanism that constantly re-centres the interface. As the temperature drops 
and more and larger bubbles of the "wrong" phase form within the main domains, 
the interface will sometimes combine with them and so we might expect it to move 
faster. In fig. 14, we do indeed see an increase in speed as the temperature drops and 
the speed appears to diverge as /? — *• /3c- This is shown more clearly in fig. 15 where 
we fit the speed to a function of the form 

— (18) 

Our results suggest values of 

/?c = 8.14(3), f5Rw = -0.71(6). (19) 

This value for j3c is consistent with that obtained by other methods, and the dy- 
namical exponent governing the critical acceleration of the random walk is given by 



Studies of discrete models and continuum growth equations |^ in condensed 
matter systems make predictions for the behaviour of the width of an interface, 
usually defined to be the square-root of our second moment. 



W{Ly,t) = ^<i)^-(t) , (20) 

as a function of the time, t, since its formation with initial width zero. For times 
much smaller than some critical time, tx{Ly), determined by the interface size, an 
exponent (3g will govern the growth of the width 

W{t)^t''<=, t<^tx{Ly). (21) 

For times much larger than this the width is governed by 0:2, the roughness exponent, 
as we discussed in section 5. Fig. 16 shows an example of the growth of a 2'(3) 
interface. At high t the size saturates as usual, but at low t we see growth governed 
by an exponent, as in the condensed matter models. For our interface, the average 
profile yields 

= 0.61(2). (22) 
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Interface Growth (Example) 
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Figure 16: An example of the growth in size of an interface for (3 = 8.5, Ly = 36, 
and a growth profile averaged over all Ly and (3. In each case, (3g is estimated from 
the slope of the dotted line. 
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6 Conclusions 



Despite the increasing fluctuations and decreasing height of the interface as the 
critical temperature is approached we have shown that it is quite practicable to 
isolate it and to measure its geometrical properties. As one might expect of a second 
order transition the fluctuation moments diverge as /3 | /3c in a manner which is to 
a good approximation independent of the transverse size of the lattice, Ly. That 
the relation between the three moments measured is more or less as predicted by 
the simple model (|^) may be fortuitous because the model does not describe the 
dependence of the moments on the transverse size very well. One of the reasons 
why we are able to measure the fluctuations reliably is that the intrinsic width 
remains small; from the measurements we have made there is no real indication 
that this quantity does not remain finite as [3 [ Pc- The dynamical (in computer 
time) behaviour of the interface is consistent with the usual behaviour encountered 
in other systems. 

From this work and [0, ^ it seems clear that the Z{N) interface in the Euclidean 
formulation of the theory really does have all the properties expected of such an 
interface in an ordinary statistical mechanical system as well as being consistent 
with weak coupling perturbation theory. The only objection to such a conclusion 
is that most of this work has been done at Lt = 2. However the behaviour of 
the interface free energy in SU{2) at high temperatures has been studied in for 
Lt = 2, 3, 4, and 6 and they find essentially no Lt dependence. 

The relation of this picture to Minkowski space is as yet unclear. In the Euclidean 
path integral, the different Z{N) vacua, and consequently domain walls between 
them, will appear in the ensemble of gauge field configurations that is integrated over 
to yield expectation values. We expect at least some observable quantities, call them 
Q, to depend upon the presence of these vacua in the sense that if there were only one 
vacuum Q would be different. The observables must be the same when calculated 
in Minkowski space as in Euclidean space. It follows that there must be some 
(necessarily non-perturbative) feature of the ensemble of gauge field configurations 
in Minkowski space which is responsible for Q taking a different value from the 
single Euclidean vacuum prediction (or else there are no Qs). The Q observables 
potentially correspond to interesting physical effects such as baryogenesis and CP 



violation in the early universe ||25|, these calculations can be criticised on the 
grounds that they assume the Euclidean domain structure carries over to Minkowski 
space but it may be that the calculation can be reformulated entirely in Euclidean 
space in which case the results should be unambiguous. 
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